Quasiperiodic Atom Optics, Focusing, and Wave packet Rephasing 
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We propose a laser field configuration which acts as a quasiperiodic atom optical diffraction 
grating. Analytical and computational results for the atomic center-of-mass wavefunction after the 
grating reveal a quasiperiodic density pattern, a semiclassical focusing effect, and a quasiperiodic 
self-imaging of the atomic wavefunction analogous to a Talbot effect. 



Increased attention has been given to periodic atom optical systems and to drawing parallels between solid state 
physics and the corresponding atom optics in an easily-controlled environment |L| . In the transient interaction regime, 
C<~) \ periodic atom optical elements have been used for spatial- 0] and time-domain ^|J4| interferometry, atom focusing and 
lithography || , and the Fresnel self-imaging (Talbot effect [g]f) of an atomic wave packet . In addition, optical lattices 
& have been used to simulate solid state effects using atomic de Broglie waves. These atoms exhibit quantized motion 
|)|, extended/localized state transitions Jl(j, and Bloch oscillations pj| . To further the connections with condensed 
O |' matter and classical optics, the extension of atom optical experiments to quasiperiodic systems seems natural [p^|-^4|. 

Recently, Guidoni and coworkers Jll| trapped cesium atoms into a quasiperiodic optical lattice formed by a three- 
dimensional laser configuration with incommensurate spatial components of intensity. Qian Niu and coworkers have 
been analyzing such optical lattices in one- and two-dimensions to understand their eigenstructure We propose 
to extend atom lithography and interferometry experiments to create atomic beams or cold trapped atoms with 
quasiperiodic center-of-mass wave functions. We present analytical and computational results for the atomic wave 
t^j , function and density after interaction with a one-dimensional, quasiperiodic atom optical diffraction grating. Mo- 
mentum distributions reveal the quasiperiodic nature of the wave packets. Furthermore, we show the possibility for 
wave packet revivals, essentially a quasiperiodic Talbot effect, where the initial atomic wave function can be (nearly) 
• recovered. 

Oh| A schematic of an experiment, similar to typical atomic focusing and Talbot arrangements, is shown in Fig. 1 
||[7]]. Two pairs of off-resonant laser beams of width w intersect at a point along a monoenergetic, transversely-cooled 
,_H ' atomic beam with velocity v z propagating in the z-direction (or alternatively in the y-direction). Assuming the laser 
. beam pairs are detuned from one another, we can ignore the cross-terms in the intensity pattern. The atoms are 
modeled as two- level center-of-mass plane waves with upper-state lifetimes of T . Hence, this optical configuration 
^sl" , forms a light-shift potential for the atomic ground state, 



V(x, z) = Vi (z) cos 2kx + V-2 (z) cos V2kx, (1) 
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where V n (z) ~ H |J7 n (z)| /8A„ for Rabi frequencies Q n (z) and atom-field detunings A„ of the laser beam pairs n = 1, 2. 
Spontaneous emission during the atom-field interaction is ignored by assuming that |f2„(z)| 2 T?«/4w z A^ <C 1. The 
quasiperiodicity arises from the incommensurate wave vectors of the optical potential, and phase stability between 
" ^ • the laser fields is essential for the potential's integrity. 

The atomic motion is described by an effective Schrodinger equation for the transverse wave function (f>(x, t) in the 
■ atomic rest frame, z = v z t. Assuming l/2Mv 2 3> (\/2Mv 2 + l/2Mv 2 + V{x, z)), this equation reads [ p7[ 
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The potential in the interaction region appears as a pulse of duration r = w/v z in the atomic rest frame. The pulse 
" ' shape is determined by the transverse laser profiles. 

The discussion here is restricted to the Raman-Nath regime for square pulses, V(x, t) — V(x) for — r < t < and 
V(x, t) = for all other times, where 

V(x) = Vi cos 2kx + V 2 cos V2kx. (3) 

In the Raman-Nath approximation one assumes that u>k 1^(^)1 t 2 /h <C 1, where u>k — 2hk 2 /M is a two-photon recoil 
frequency; the kinetic energy is ignored during the interaction, allowing an immediate integration of Eq. (^|) using 
(f>(x, -t) = 1 and Eq. (§) Q. We can write <j)(x,0) = exp [-iV(x)r/h] or 
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</)(x, 0) = exp i(A\ cos 2kx + cos y2ka 



(4) 



where the pulse area A n = —V n T/h. The standing-wave light fields acts as a quasiperiodic atomic phase grating. 

To follow the evolution of the wave packet after the interaction, we can expand the exponentials in Eq. in 
a plane wave representation. The resulting, initial wave function for the time-dependent Schrodinger equation is a 
superposition of free particle eigenstates exp[ip mtn x/h] with energies E m ^ n = pl^ n /2M, giving the result 



(x,t> 0) 
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is the phase of the momentum component 
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Pm,n = 2Hk(m + n/V2) 
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and J m is a Bessel function of the first kind. The momentum space wave function superposes a set of regularly 
spaced, but not periodic, components which are integer combinations of momentum exchanges between the atom 
and fields. By squaring Eq. (^|) and using a sum rule for Bessel function products, the transverse atomic density, 
p(x,t) — <p*(x,t)4>(x,t), can be written as 



p(x, t) = 22 p m ,n(ukt) exp 
The density has time-dependent Fourier amplitudes 
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Pm,n{ukt) = J m \ sin[(m + —)LU k t] 



x J„ I 2A 2 sin[(— = + -)uj k t] 



V2 2' 
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creating a spatial pattern which evolves in time. Thus, the atomic density is not only a quasiperiodic function of the 
coordinate x, but of the coordinate z = v z t as well. 
The density Fourier transform (DFT), 



P(Q,t) 



— p(x, t)e~ %qx = ^ Pm,n(ukt)8(g - Pm,n/fr), 



(10) 



has peaks at q = 2k(m + n/y/2) by Eq. (^). When squared, p(q,t) gives the time-dependent structure factor of 
the atomic distribution. In realistic experiments the delta-function lineshape of each spectral component would be 
broadened according to the initial momentum distribution of the transverse atomic beam. For a thermal velocity 
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(uj k t) exp - (p m , n ut/2hf in Eq. @ 



distribution with most probable speed u, the replacement p„ 
is sufficient to account for Doppler dephasing. 

To detect the density as a function of t, one can scatter a transient probe off of the atoms at t to record the 
time evolution of certain Fourier components of the density. For example, a weak probe pulse with duration < 

-kx backscatters a field Ej, s proportional to p\fi(ui k t) in the +a?-direction: E^ s ~ 



(ku) and wave vector k p = 



(hut) . This is a type of free induction decay experiment to detect 



Ji (2Ai sin[uj k t}) Jo (2^2 sm[u>kt/y/2]) exp 

ground state population gratings jis|| . 

More importantly, either direct atomic deposition or lithography using the atomic beam to impinge on a prepared 
substrate would reconstruct the atomic density at a fixed time. Atomic lithography has advanced to the point where 
atoms can carve nanostructures in materials such as silicon, silicon dioxide, and gold pC( |. Such quasiperiodic surfaces 
could be used for solid state surface and transport studies. The implications for quantum and optical properties, 
including photon localization, may be profound owing to the quasiperiodic boundary conditions for the electron or 
optical waves Q . 
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We now examine two different phenomena in the pattern formed by the atoms (§1), ignoring dephasing. Semiclassical 
(near-field) dynamics explain a focusing effect, similar to that seen after periodic phase gratings f|. In Fig. 2 the 
optical potential is shown for A\ = 5 and Ai = 10. Each potential well acts as a lens which can focus atoms 
using the impulsive (dipole) force, Fix) — MAv(x)/t = —dV(x)/dx, where Av(x) is an impulsive velocity kick. To 
illustrate this effect, V(x) is Taylor expanded around its minimum at x — to give the focusing force near this point, 
F(x) = MAv(x)/t w (4Vi + 2V<2)kx. Solving for Av(x) and setting tf — x/Av(x), this geometrical argument yields 
a focus at the time u>ktf = (2A± + A2)~ 1 that translates into a spatial distance Zf — v z tf. The ratio of pulse areas in 
Fig. 2, A2/A1 — 2, was chosen so that each standing wave contributes an equal semiclassical force. 

The density at this " quasi" -focus with its peak at x — is also shown in Fig. 2. Additional density peaks result 
from focusing by the shallower wells which occur at the quasiperiods of the potential. For example, the peak at 
2kx w 7 * 2tt ps 5y/2 * 2n occurs near a potential well where the incommensurate standing waves are nearly in phase. 
In general, the irrational wave vector ratio, \[2 in this case, can be approximated as the ratio a s /b s , where a s and b s 
are positive integers without common factors. Quasiperiods will then be defined by 2kx rs 2irja s ~ 2'Kjb s \f2 for any 
integer j. A converging sequence G s which approximates \[2 is given in Table 1 p2[. 
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2 


3 


4 


5 


6 


a s 


3 


7 


17 


41 


99 


b s 


2 


5 


12 


29 


70 


G s = a s /b s 


1.5 


1.4 


1.4167 


1.4138 


1.4143 


V2-G s 


8.6e-2 


1.4e-2 


2.5e-3 


4.2e-4 


7.2e-5 


u k t s 


47T 


20tt 


24tt 


116tt 


1407T 



Table 1. Sequence of approximations for \[2 defines the 
spatial quasiperiods, 2kx w 2nja s , and the quasi- Talbot 
times t s in Eq. (O). 



Several peaks in Fig. 2 are labeled by their values of (ja s ,jb s ). For smaller differences between ja s and jb s y2, 
the quasiperiodicity is more pronounced (i.e., the density peaks near 2kx « 2tt ja s approach the size of the peak at 
x = 0). Thus, the density reflects the quasiperiodicity of the system. 

The DFT contains the spectral information important for lithography or scattering experiments. In Fig. 3 we plot 
the Fourier amplitudes \p m ,n{uJkt)\ [Eq. (0)] at the wave vectors q = 2k(m + nj \/2) > and at t = tf for A\ = 5 and 
A2 = 10 again: this is the DFT of Fig. 2. The major peaks are labeled by (to, n) to show the variation of amplitudes. 
The inset of Fig. 3 has the same axes as the main graph and shows the range of q values with significant amplitudes. 
This Fourier spectrum has qualitative scaling properties: for any Fourier wave vector of the density qo, a wave vector 
q' = 2fc(m' + n'/V2) can be found which is arbitrarily close to go, even if the amplitude of that component is much 
less than one. 

Atoms that propagate after interacting with periodic atom optical elements exhibit a self-imaging of their wave 
function - Talbot effect - and return to uniform density at times w^t — 2irj for integer j > While an exact 

self-imaging which reproduces 4>(x,0) [Eq. ([5])] is impossible owing to the dispersion in Eq. (6), the quasiperiodic 
nature of the wave function can lead to a quasi-rephasing when ip m ^ n ~ 2irj' for some integer j' > 0. 

The appropriate observation times for self-imaging will produce phases for each momentum component which are 
nearly integer multiples of 2n. For odd values of b s , the choice, uJkt s = 47r6 s , from Eq. (0) gives the phases 



<Pr, 



l (47r& s ) = 2ir(2b s m 2 + b s n 2 r, 



(11) 



where we again refer to Table 1, and b s \/2 ~ a s by construction. For even values of 6 S , the rephasing occurs at 
UJkts = 2irb s from Eq. ([Tl|). At these times the first two terms in Eq. ( |i"T|) are integers for all m, n. Furthermore, the 
third term is nearly an integer, as required. 

The wave function phase 9(x,t) is defined by <f>(x,t) — \cf)(x, t)\ exp[i6>(x, t)]. For exact self-imaging, the wave 
function should be the unitary exponential, Eq. (Q), which has density equal to one and phase 9(x, 0) = A\ cos(2fcir) + 
A2 cos(V2kx). In Fig. 4 we plot 6{x,t B ) for A\ = A% = 1 and s = 2,4, and 6. The self-imaging becomes more 
pronounced at longer times t s , corresponding to a better approximation of \[2 by G s . The average values (denoted 
by the bar) and standard deviations (denoted by a) of both p(x, t s ) and the phase difference, 5 S = 8(x, t s ) — 6(x, 0), 
are shown in Table 2 for the cases of Fig. 4. 
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s 


2 


4 


6 


P(x,t s ) 


1.019 


1.00085 


1.000024 


a(p(x,t s )) 


0.505 


0.0170 


0.000508 


S s (rads) 


0.0324 


0.00645 


0.00106 


a(S s ) (rads) 


0.417 


0.0899 


0.0156 



Table 2. Average and standard deviation of 

p(x, t s ) and 5 S = 6(x, t s ) - 6(x, 0) for \2kx\ < 16tt. 

The improvements in the density and phase are evident as a(p{x,t s )) and cr(5 s ) converge monotonically to zero as t s 
increases. Furthermore, a focusing effect, like the one detailed above, will occur near the times t s + 1/ since the wave 
function at t s is nearly identical to (f>(x, 0). 

Of course, the ability to perform an experiment on long time scales is severely limited by transverse cooling 
considerations. An initial momentum width of Ap x in the atomic beam washes out these coherent effects in a time 
t < XM / Ap x ~ (fcu) -1 , the Doppler dephasing time. This condition does not present a problem for thick lens focusing 
and lithography schemes |5|,^(|. In order to see these thin lens effects, the atom beam must be cooled or collimated 
near the recoil limit for focusing (ku < LOk) or below the recoil limit for Talbot self-imaging (ku <C Wfe). Recent 
experiments have conformed with the focusing condition PUp^ . 

In summary, this letter has introduced the possibility of quasiperiodic atom optical elements made from laser inten- 
sity gratings with incommensurate wave vectors. The analytical results show that the atomic wave packet becomes a 
quasiperiodic function, developing momentum components which are similarly incommensurate and therefore regular, 
but not periodic. The atomic density is a function of the time of flight from the diffraction grating. Atoms come to 
semiclassical "quasi" -focuses according to the depth and curvature of the potential wells. The quasiperiodic density 
pattern can be used to create a quasiperiodic surface for condensed matter studies when used for atomic lithography. 
Furthermore, ultracold atoms will exhibit a quasi-self-imaging of the wave function. 
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FIGURE CAPTIONS 

Fig. 1. The atomic beam traverses the quasiperiodic potential formed by the laser beams and is detected by light 
scattering or lithography after free propagation. 

Fig. 2. The optical potential in the interaction region for A± = 5 and A 2 = 10 and the corresponding atomic 
density at the quasi-focus, Zf — v z [oJk(2Ai + A2)] . Quasiperiods are pronounced near 2kx ~ 2wja s « 2wjb s \/2, 
where peaks are labeled by (ja s ,jb s ). 

Fig. 3. Fourier amplitudes \p m , n ( u >k'tf)\ at q — 2k(m + n/\^2) for the density in Fig. 2, as given by Eq. (||). Larger 
amplitudes are labeled by (m, n). Inset: Full range of Fourier components with significant amplitudes. 

Fig. 4. Quasiperiodic Talbot Effect for A\ = 1 and A2 = 1. We plot 9(x,t s ), the phase of the atoms at the quasi- 

Talbot times t s for s = 2 (-.-.-), 4 (....), 6 ( ), versus 8(x,0), the initial phase ( — ), where 9(x,0) — Ai cos(2fcx) + 

Ai cos 
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